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HYDRODYNAMICAL SPECTRAL EVOLUTION FOR RANDOM 

MATRICES 


PETER J. FORRESTER AND JACEK GRELA 


Abstract. The eigenvalues of the matrix structure X + where X is a 

random Gaussian Hermitian matrix and X^*^^ is non-random or random inde¬ 


pendent of X, are closely related to Dyson Brownian motion. Previous works 
have shown how an infinite hierarchy of equations satisfied by the dynamical 
correlations become triangular in the infinite density limit, and give rise to 
the complex Burgers equation for the Green’s function of the corresponding 
one-point density function. We show how this and analogous partial differ¬ 
ential equations, for chiral, circular and Jacobi versions of Dyson Brownian 
motion follow from a macroscopic hydrodynamical description involving the 
current density and continuity equation. The method of characteristics gives 
a systematic approach to solving the PDEs, and in the chiral case we show 
how this efficiently reclaims the characterisation of the global eigenvalue den¬ 
sity for non-central Wishart matrices due to Dozier and Silverstein. Collective 
variables provide another approach to deriving the complex Burgers equation 
in the Gaussian case, and we show that this approach applies equally as well 
to chiral matrices. We relate both the Gaussian and chiral cases to the asymp¬ 
totics of matrix integrals. 


1. Introduction 


One of the most basic questions in random matrix theory asks for the limiting 
global spectral density, given the distribution on the elements of the matrices or the 
distribution on the space of matrices. Perhaps the most celebrated result in this 
class is the Wigner semi-circle law. It applies to real random symmetric, or complex 
Hermitian matrices, in which the diagonal entries are independently distributed 
with mean zero and variance unity, and the upper triangular off-diagonal entries 
are independently distributed with mean zero and variance two. Random matrices 
of this type are referred to as Wigner matrices. Consider a specific class of Wigner 
matrices of size N x N. A scaling of the matrices so that the eigenvalue density has 
compact support is referred to as a global limit. For Wigner matrices this is achieved 
by dividing each matrix by \/2N before taking the A —>■ oo limit. Moreover, the 
corresponding global spectral density normalised to integrate to unity, is 

given by 



|x| < 1 
|a:| > 1. 


(l.I) 


In its fully generality, the most common proof of this result, and in fact the one 
provided by Wigner [33] , is to show that 



( 1 . 2 ) 
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where the RHS is the fc-th Catalan number. This is equivalent to establishing that 
the moments jdx are given by the fc-th Catalan number for each 
k = 0,1,2,.... The functional form (HI) now follows as a consequence of the 
Wigner semi-circle law being the unique probability density with such moments. 

An extension of the question seeking the global density of a single Wigner matrix 
is to ask for the global density of the sum of matrices 

Y = X^°'>+X, (1.3) 

where X is again a Wigner random matrix, but may be non-random, or random 
independent of AT. To answer this a tool kit beyond the analysis of moments is 
required. Specifically, the analytic properties of the Stieltjes transform plays an 
essential role. Thus let X'^^'>/\/2N have global spectral density and define 

the corresponding Green’s function (or Stieltjes transform) by 

G«»(,)=/ (G4) 

Jl(0) Z-X 

where is the support of Then, with G{z) denoting the Green’s function 
of the sought global density p{y) of Yj\/2N, it is a known result (see e.g. [151 
Th. 18.3.2]) that G{z) is determined as the solution of the functional equation 

G(z) = GW(z-iG(z)), (1.5) 

subject to the requirements that G{z) ^l/zas| 2:|—>-00 and that G{z) be analytic 
for z ^ I, where I is the support of p{y). 

As an illustration of m ,suppose 

p(0)(x) = (5(x), (1.6) 


the Dirac delta function at the origin, in which case G^°^(z) 
functional equation (H reads 


G{z) 


1 

z-G(2)/4’ 


and this has the solution 


1/z. Then the 


G(z) = 2{z — \/ z^ — 1). 


(1.7) 


But according to the Sokhotski-Plemelj formula for the inverse of the Stieltjes 
transform, 

and substituting dnD reclaims dnj. 

Our interest is to develop a viewpoint of the functional equation (EH) for the 
Green’s function of the density for the matrix structure (ESI) as having origins in 
the hydrodynamical equation 

f^ (u (X) + rts; f) 5) <i„) ), (1.9) 

where Vi{x), V 2 {x,y) are particular one and two body potentials, and f is a scaled 
parameter. The scaling is required to compensate for the normalisation of p{y, r) 
being such that its total integral is unity rather than N. As to be reviewed in Section 
H in the case that X in (E31) is a standard Gaussian matrix, it is well known that 
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a hydrodynamical equation, in particular the complex Burger’s equation, relates to 

m- 

The advantage of viewing m as a consequence of (HI) is that a hydrodynamical 
description applies equally well to the case that X and have the chiral structure 


X = 


Omxm 

Yt 

.nxm 


Zm'Kn 

Onxn 


n > m, 


( 1 . 10 ) 


and similarly with Z replaced by where Z is a standard real or complex 
Gaussian matrix, and with Z^^'^ non-random or random independent of Z. The class 
of random matrices (|1.10l) is fundamental to random matrix theory (see e..g. [171 
§3.1]). Although a hydrodynamical description of the spectral evolution of complex 
chiral Gaussian matrices has been given in the recent works Ell, the general solution 
of the resulting partial differential equation for the Green’s function seems to have 
not been considered. In Section [3] we show how (11.91) leads to a partial differential 
equation, transformable to an inhomogeneous complex Burger’s equation. The 
partial differential equation was first derived in [H], H, with the starting point 
in the latter being a microscopic description involving Dyson Brownian motion as 
reviewed in the introduction to Section [S] Moreover, we show how to solve this 
partial differential equation in a form analogous to dlH). In so doing we are able 
to reclaim a functional equation for the Green’s function of so-called non-central 
Wishart matrices given by Dozier and Silverstein [12] . 

Specifically, consider random matrices of the form 


W= \[^Z + Z^°\ 

V m 

where Z, Z^^'> are rectangular mxn matrices, Z a standard Gaussian and Z^^'^ fixed 
or random independent of Z. Let a = lim„_m->oo(iT-/w) — 1, and let the m,n —>■ oo 
limiting density of eigenvalues of be equal to and the limiting 

density of eigenvalues of W^W be equal to p^{y,f). In this setting, and with 



1+ C i?+ denoting the support of p^{y,T), 


VVV-. vv 




g^{z\T) = Fw 


m^y) 


j( 0 ) z{Fwy - 2afTw - y 


(l.II) 


( 1 . 12 ) 


where Fw = 1 — 2Tg^ { z ; t ). With the identifications t = cr^c/2, d = (1 — c)/c and 
= —m this is Eq. (1.1) in [12] . 

When both X and X^'’^ in are unitary, we are dealing with circular en¬ 

sembles for which a hydrodynamical equation is well known [3D]. We review 
the necessary working in Section |3| and proceed to present a novel hydrodynamical 
description of Jacobi ensemble in trigonometric variables, which can be interpreted 
as a ’’chiral” deformation of circular ensembles. In Section [5| and IH we turn to the 
collective variables approach which is used to review the well-known asymptotic ex¬ 
pansion of Harish-Chandra/ Itzykson-Zuber integral and to study not previously 
considered expansion of Berezin-Karpelevich type integrals. 
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2. Gaussian ensembles 

In this section we review how the hydrodynamical equation (HH) comes about 
from the study of the eigenvalues of the matrix sum in the case that X is 

a standard Gaussian matrix with real (/3 = 1) or complex (/3 = 2) entries |13) . 
and furthermore leads to the functional equation (lESD. The starting point is the 
fact that the Gaussian distribution on the space of matrices {T} which are real 
symmetric (/3 = 1) or complex Hermitian (/3 = 2), 

F) = ^ exp ( - |-Tr(F - , (2.1) 


where denotes the normalisation constant, satisfies the diffusion equation 


dPr 

~d7 


1 




d'^Pr 


( 2 . 2 ) 


In (12.21) the label /i ranges over the independent elements, including both the real 
and imaginary parts of the off-diagonal elements if they are complex, and = 1 
for the diagonal elements, and = i for the off-diagonal elements. 

There is also a class of Hermitian matrices with distribution (EB that satisfy 
(12.211 with ,5 = 4. Thus the Hermitian matrix F, and similarly is now a 

2N X 2N matrix formed from an N x N matrix with each element a 2 x 2 matrix 
of the form 


z —w 

—w z 


z,w G C. 


Such 2x2 matrices are isomorphic to the real quaternion division algebra, one of 
only three associative real normed division algebras along with the real and complex 
numbers, so in this case F is said to have real quaternion entries; see e.g. m §i-3.2]. 

The differential operator on the RHS of (12.21) can be interpreted as the Laplace- 
Beltrami operator associated with the metric form for the matrix spaces 


{dsf = Tr(dFdFt) = ^ 5 ^,dF^dF„ 

with /i labelling the independent elements as in (1^ and similarly i/, and where 
see e.g. [HI eq. (11.9)]. Introducing the diagonalisation formula for 
F, F = ULU\ where U is the matrix of eigenvectors and L is the diagonal matrix 
of the eigenvalues, allows to be rewritten according to the separated form 





J=l[\Xk-Xjf, 

j<k 


(2.3) 


where the operator Ojj involves derivatives with respect to variables relating to 
the eigenvectors only. The significance of this is that the eigenvalue distribu¬ 
tion Pt{Xi, ..., Xn) obtained by integrating over the angles U and the distribution 


Po{X^°y) of in (lO) . 


p^(Ai,..., Aat) = J J dU J dX^^'^Pr C/Lt/I') Pq (2.4) 
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satisfies the Smoluchowski-Fokker-Planck equation 


dpr _ . 

dr 


N 


i=i 


d /dW 
dXi 



(2.5) 


with 

W = - log|A,-Afc|. (2.6) 

l<j<k<N 

In the case that the r dependence in is modified so that Pr satisfies not 
( 1 ^ . but the heat equation for Brownian motion in a harmonic potential i.e. the 
Smoluchowski-Fokker-Planck equation (12.51) with W correspondingly modified by 
the addition of an harmonic potential | first derived by Dyson [H] . 

As such the corresponding process is referred to as Dyson Brownian motion. 

As pointed out in El: the Smoluchowski-Fokker-Planck equation (12.51) has a 
more standard interpretation than its origin in random matrix theory. Specifically, 
consider a classical system of JV particles interacting on a line with potential W. 
Suppose the particles execute overdamped Brownian motion in a fictitious back¬ 
ground fluid with friction coefficient 7 , and furthermore the system is at inverse 
temperature /3. It is a basic fact — see e.g. [32] — that in this setting the evo¬ 
lution of the probability density Pt-(Ai, ..., Aat) for the location of the particles 
at positions Ai,... An is given by the Smoluchowski-Fokker-Planck equation (12.51) . 
where the LHS is to be multiplied by the friction coefficient 7 . The random ma¬ 
trix problem gives rise to the specific potential ( 12 . 6 p . corresponding to the particles 
interacting via the repulsive pair potential V 2 {x,y) = — logja: — y\, and so the 
underlying classical gas is referred to as a log-gas m 

Our interest is in the one-body dynamical density p{x; r) defined as an average 
over (I37T]) . 


=-^ l'^S{x-\i)\ , (2.7) 

\i=l / 

which has been normalized to have total integral unity. We probe a global regime 
for which the eigenvalues are scaled so that they have finite support as in dm. 
Note that since the integral over x of p{x; r) is unity, the integral of Np(x; t) over 
X is N. In this scaling the inter-particle spacing goes to zero, and the response of 
the system to perturbation is governed by macroscopic equations m- Relevant to 
the one-body dynamical density is the macroscopic equation 

N^J{x-,t) = T{x]t), (2.8) 


where N^J{x-,t) is the one-body current related to the density by the continuity 
equation 

■^p{x]f) =--^J{x]f), T = Nt (2.9) 

(the scaled parameter f is introduced to compensate for the integral of p{x', r) 
normalised to unity) while P{x]t) refers to the macroscopic force density. For the 
log-gas in the long wavelength regime the force density to leading order is of an 
electrostatics origin, so implying (see e.g. [T3] l 


d ( f 


'OO 
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Differentiating both sides with respect to a;, and making use of the continuity 
equation (12.9p on the RHS, we see the hydrodynamical equation (11.91) results with 
Vi{x) = 0 and V 2 {x,y) = -log|a; - y\. 

We now want to show how this particular hydrodynamical equation leads to the 
functional equation m- For this purpose, it is convenient to introduce the Hilbert 
(or Cauchy) transform as the principal value integral 

'H[n](x):=PV [ dy, xGl. (2-11) 

Ji x-y 

The hydrodynamical equation of interest then reads 

We also introduce the Green’s function 

G(z;f) := f dy, (2.13) 

Ji z-y 

cf. (11.41) . 

Next we follow the working in §III.B.4], which begins by noting that as a 
consequence of the residue theorem, the Green’s function is related to the Hilbert 
transform by 

G±(x;f) = =F*7rp(x;T)+'H[p(-;f)](x), x G/. (2.14) 

Using this in (12.121) gives 

2^(G-(a;;T) - G+(x;t)) = -^((G_(x;f))2 - (G+(x;f))2). (2.15) 

It must therefore be that the function 

2^G(z;f) + ^(G(z;t))2 

is analytic throughout the entire complex plane. But according to (I2.13L G(z; f) ^ 
1/z as \z\ —>■ oo, so this function furthermore goes to zero at infinity. The only 
analytic function with this property is the zero function, and so after minor ma¬ 
nipulation we have 

o o 

—G(z;f)-b G(z;f)—G(z;f) = 0. (2.16) 

OT oz 

This is the Euler equation, or equivalently complex Burgers equation of hydrody¬ 
namics. Thus with z = X + iy and G{z-, t ) = U + iV, {U, V) is the velocity field at 
point (x, y) in the plane for an ideal fluid at constant pressure. 

To solve the initial value problem for (12.161) . we invoke the method of complex 
characteristics where both z and G are complex functions. This is a slight gener¬ 
alization of a standard technique applicable to inital value problems of real first 
order PDEs. We present this method for a general first order equation of the form 

A(G{z;f),z,T'^-^G{z;f) + b(g{z;t),z,t^-^G{z;t) = g(^G(z; f), z,f). 

(2.17) 

The main idea is to seek a coordinate transform (z, f) —>■ (a, /3) such that the PDE 
(12.171) becomes an ODE along the curves of constant a, 

^G(z; f) = g(g(z; f), z, f) 


( 2 . 18 ) 
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called characteristic lines or simply characteristics. By the chain rule ^ + 

the left-hand sides of (I2.18P and (12.171) dictate the system of equations de¬ 
scribing the characteristics, 

^f(a, P) = A{G{a, /3), z(a, /?), f(a, , 

^z(a,/3) = (2.19) 

where G(a,/3) = G(z(a,/3);f(a,/?)). These equations form lines in the (z,f) 
space, labeled by the fi parameter and passing through the prescribed initial point 
(z(a, 0), f(a, 0)). For the latter to be determined, the Green’s function G{z -f = 0) 
on the f = 0 line is required. With this initial data specified, the set of equations 
(12.181) and (12.191) are in principle solvable by standard means and comprise the 
sought solution to (12.171) . 

The equation (I2.16P is an instance of (I2.17P with A = 1, B = G and G = 
0. We read off from ()2.17p and ()2.19p that the differential equations describing 
characteristic lines and the propagation of the solution are 

^z(a,/3) = G(a,/3), ^f(a,/3) = l, ^G(a,/3) = 0. (2.20) 

The initial data comprises of the initial position z(a,0) = a,f(a,0) = 0 and the 
starting Green’s function G(a, 0) = G( 2 (a, 0); f (a, 0)) = G(a;0). Explicit integra¬ 
tion gives 


G(a,/3) =G(a;0), 

f(a,/3) = /3, 

z(a,/3) = a-b/3G(a;0). (2.21) 

These, after eliminating a and /3, yield the functional equation 

G{z;t) = G(^z-TG{z;T)-0y (2.22) 


which is to be compared to (11.51) . Equivalently, recalling the definition dH , we 
obtain the implicit integral equation 


Jj( 0 ) z - tG(z;t) - fj, 


(2.23) 


Working closely related to the above discussion can be found in [33] . 

To anticipate the precise relationship between (12.221) and (11.51) . let us choose 
T = 1/(4A^). The RHS of (12.11) is then proportional to exp{—NP{Y — and 

thus we see that Y = X + where X = X/\/2N with X a standard Gaussian 
matrix. This is precisely the setting which gives rise to (IE3. On the other hand, 
the choice r = 1/(47V) is, according to (12.91) . equivalent to the choice f = 1/4, and 
this substituted in (12.221) gives (11.51) . 

In the Introduction the functional equation (II.5|) was illustrated by showing that 
the case p^°\x) = 6{x) leads to the Wigner semi-circle law (11.11) . Another example 
which permits an explicit functional form for the density is when 


G(°)(z) = i( 


1 


z — a z Y a 


(2.24) 
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This corresponds to an initial density 

- a) + 5(x + a)), (2.25) 

or equivalently to in (lESJ being a diagonal matrix with half its eigenvalues at 
a and the other half at —a. Substituting (12.2411 in (12.221) with f = 1 shows that the 
Green’s function G{z) satisfies the cubic equation 

- 2zG^ + {1 - + z'^)G - z = 0. (2.26) 

With ^ = z — G, this equation first appeared in the present context in [S] , where it 
was shown to correspond to a spectral density supported on two disjoint intervals 
symmetrical about the origin for o > 1. In the case a = 1 the intervals meet at the 
origin, and it can then be shown that the eigenvalue density has the explicit form 
[291 eq. (6.118)] 

p{y) = ((3^3 + v'27-8y2)2/3 _ (3^ _ (2.27) 

for —y^27/8 < y < jS. For a discussion of (12.261) in terms of caustics corre¬ 
sponding to the complex Burger’s equation, see [33] . 


3. Chiral Gaussian ensembles 


3.1. Partial differential equation. We now turn our attention to the matrix 
sum mi in the case that X is a chiral Gaussian random matrix as specified by 
(I1.10|l . and has the same block structure as X but is non-random or random 
independent of X. Thus the matrix Y in (11.31) similarly has the block structure 


Y = 


'-'mxm 
. nxm 


IWnxn 

Onxn 


W = Z + Zo 


(3.1) 


Instead of (EtD, we now have a distribution on the block matrix W specified by 
Pr{Z^°'>;W) = -^exp(- -^Tr(lF-Z(°))t(lF-Z(°))). (3.2) 

The diffusion equation but with Y replaced by W, again allows for a 

characterisation of this matrix distribution. Before the (Hermitian) matrix Y was 
decomposed according to the diagonalisation formula. The appropriate decompo¬ 
sition of the matrix W is the singular value decomposition W = ULV\ where U 
and V are real orthogonal (/3 = 1) or complex unitary {j3 = 2) matrices of size 
m X m and n x n respectively, while L = diag(a;i,..., Xm), where {xj} are the sin¬ 
gular values of W or equivalently {x|} are the eigenvalues of W^W. As discussed 
in [MU], integrating over the distribution of Z^°'> leaves a distribution function 
Pt{xI, ... ,x^) depending on the parameter r and the eigenvalues W^W only. To 
characterise this distribution as the solution of an evolution equation, we require 
the fact (see e.g. [T71 §11.2.2]) that the Jacobian J in the formula (12.31) should now 
read 

m 

n a = n-m+l-2/l5, 

j—1 l<j<k<m 

and furthermore on the RHS the replacements 


{Aj} I—>■ {xj}, N ^ m 


( 3 . 3 ) 
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should be made. Doing this then gives that Pt{x\, ..., satisfies the Smoluchowski- 
Fokker-Planck equation (12.51) with 


/ m 

W = \og\xl-x% (3.4) 

j—1 l<j<k<ra 

where a' = a + 1//3 and with the replacements (13.31) . In the log-gas analogy, the 
domain is now the half line x > 0, and there is both a one and two body potential 
given by 

Viix) =-'^^ogx'^, V 2 {x,y) =-\og\x^-y^\. (3.5) 

We now turn our attention to the hydrodynamical description of the global 
density p'^(x; r) which is defined as 


p'=(x;t) = — {Y 
m \ 

\i=i 


{5{x - Xi) -I- <5(x -b Xj)) 


(3.6) 


where we average over the measure (1^ . and the superscript “c” denotes the chiral 
case. This density is an even function in x and is normalized so that integration 
over the positive half line x > 0 gives unity. As in the discussion of Section [51 to 
access the global regime we must scale the parameter r = mr, and also scale 


a = 


lim-1, 

n^oo TO 
m—>-oo 


(3.7) 


which so determines the limiting ratio n/m. In terms of these scaled parameters, 
from the explicit form (13.5|) of the one and two body potentials, and the fact that 
the domain is a half line, the hydrodynamical equation (fm reads 


dp'^jx-^r) 

dr 


^^p‘=(x;f)^(- ^logx^ - p'=(j/; f) log |x^ - 2/^1 

^^p‘=(x;f)^(-^logx^-y p^{y]T)\og\x-y\dyj'^, (3.8) 


where the second line follows by writing log |x^ — y^| = log |x — y| -I- log |x -I- y| and 
the fact that p'^(y; t) is even in y. 

Introducing the Hilbert transform as defined in (12.111) . (13.81) can be written 

+ P^(-; r)](x)) . (3.9) 


Introducing too the Green’s function 

G\z-,t)= j ^-^dy, (3.10) 

Ji+u-u ^~y 

where 1+ C M''" is the support of p^{y,T) on the positive real axis, we see that 


G±(x;f) = =F*7rp'’(x;r)-b'H[p‘’(-;f)](x), x G/+U-/+, (3.11) 

where is defined according to (12.141) . Proceeding as in the derivation of (12.161) . 
it follows from the use of (13.111) in (13.91) that 


t) - f) + (g 9»; f) + 1 ) = 0. 


(3.12) 
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3.2. Solution of the partial differential equation. We now seek the general 
solution of the initial value problem for this partial differential equation, first in the 
case a = 0, then in the more difficult case d > 0. 

The case d = 0. In the case d = 0 we see that (13.121) reduces to the Euler equation 
(I2.16|) . Thus as with (j2.22L the solution in this case is 

G"(z;f) = G"(z-fG"(z;f);0), (3.13) 

but we must keep in mind that G'^ is defined by (13.101) rather than (12.131) . For 
(12.221) and (j3.13l) to imply identical, up to a scale factor, eigenvalue distributions 
and singular value distributions respectively for general r, we see that we must have 
the initial conditions related by 

p(2/;0) = ip^(y;0), (3.14) 

where is the support of p^{y;0) on the positive real axis. The factor of ^ is to 
compensate for the normalisation of the LHS being such that integration over the 
whole real line gives unity, while on the RHS integration of p‘^(y;0) over the half 
line y > 0 gives unity which is evident from the definition (13.61) . 

Specifically, we see that with the initial conditions related by (13.141) , the solutions 
of (12.221) and (13.131) are related by 



and thus 

p(2/; ^) = (^2/; 0 • (3-15) 

In relation to the initial condition p{y;0) = S{y), after recalling that the Wigner 
semi-circle law ()l.ll) results from the parameter value r = 1/4, we see from (13.151) 
that 

P^(x;t) = —-V^St — , 0 < a; < 2\p2^. (3.16) 

27rr 

It must therefore be that for = Omxn in (13.21) the global density of the singular 
values of the m x n, n > m, matrix W with distribution specified by (EH), is in the 
case that limm->.oo = 1 equal to this functional form. 

To see this latter result, which is well known, first note that each singular value x 
of W is related to an eigenvalue y of W^W hy x"^ = y. Changing variables according 
to this prescription in (13.161) gives the density function 

. ! yj9,T-y, 0<y<8T. (3.17) 

dTTT^ 

With T = I this specifies the Marchenko-Pastur law for the limiting density of 
the eigenvalues of the scaled matrices ^W^W, with W now a standard Gaussian 
rectangular matrix, again in the circumstance that limm-^-oo ^ = 1 ) see e.g. m 
Section 3.4.1]. This scaling is consistent with that implied by (13.11) . with = 
Omxn and r = r/m = 1 /( 8 to). 

An analogous discussion holds for the initial density (12.241) . In the case a = I, 
T = 1 / 2 , and upon changing variables y'^ = x, we conclude that the limiting density 
of the eigenvalues of the scaled matrices ;^(W -|- Z^^'>)^{X + where AT is a 

standard Gaussian m x n rectangular matrix is an m x n matrix with half its 
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entries on the diagonal equal to +1, the other half —1, all other entries equal to 0, 
is equal to 

1 (3^3 + V27 - 4x)2/3 - (3^ - V27 - 

25 / 331 / 2 ^ 2 . 1/3 ’ ^ 1 

supported on 0 < a; < This is the density function for the Raney distribution 
with parameters p = 3, r = 2 [31] . where for general 0 < r < p and p > 1 the Raney 
distribution is characterised by its k-th. moments according to 

= + fc = 0,l,2,.... (3.19) 

Note that i? 2 ,i corresponds to the Catalan numbers; recall (HU- 


The case d > 0. For nonzero parameter d, the solution of (I3.12|) to an initial 
value problem is given by the method of characteristics described in Section [21 The 
chiral ensemble is an instance of (12.171) with A = 1,B = + a/z and C = . 

Accordingly, the system of ODEs describing both the characteristic lines and the 
propagation of Green’s function is 


Az(a,/3) = G^(a,/3) + ^^, 

dp z[a,p) 

jLrAn, R\ - 

G [OL^Pj (X . ; 

dp z{a^pY 


dp 


T{a,l3) = 1 , 


with initial conditions z{a, 0) = a, f (a, 0) = 0 and G^{a, 0) = G^{z{a, 0); f{a, 0)) = 
G‘’(a;0). The last equation for t is simply solved as f = /3 whereas the first two 
are coupled but readily solved to give 


z{a,P) = Gl{a)l3\la + G^laJ/S + 2d^, 

\J a + Gl{a)l3 + 2 d^ 


f(a,/3) = /3, 


G^{a,P) = Gl{a)- 


\Ja + GQ{a)P 


(3.20) 

(3.21) 


where we used a simplified notation Go(q!) = G‘^{a, 0). Now we make the substitu¬ 
tion T = fj and calculate auxiliary formulas by multiplying the equations for z and 
G'’ and squaring the equation for z, 


G%z;t) GS(a) 

z a-|-Go(Q;)T’ 

z^a = (a -I- GQ(a)f)^a -I- 2dT(a -I- GQ(a)f). 


(3.22) 

(3.23) 


From (13.221) we find Gg = and from that 

« + (3.24) 

where from now on we suppress the arguments of G'’ and Gg for brevity. We plug 
(I3.24|) into (13.2311 to obtain a formula for , 


a 


2 


{z - G"t)2 


2aT{z — G^t) 


z 


(3.25) 
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Next we recall the definition (jd.lOl) to determine Gq and use the symmetry of the 
initial spectral density to rewrite 


GS(a)= 


AO), 


r(0) 


a — p 


= 2a 


This variant of Gq plugged into (13.221) gives 


p^’^°\p)dp 

j ( 0 ) 0-2 _ ^2 


p^’(°Hp)dp 


Jpo) 


a 


(3.26) 


(3.27) 


where we also used (13.241) . Lastly, we recall the equation (13.251) for o? and so obtain 
an implicit integral equation 


G=(z;f) = 2 


p‘^d^'>{p)dp 


z — G'^{z; t)t — 


2af 


(3.28) 


z-G^(z-f)f 


We can check, upon recalling (I3.10L that in the case a = 0 (13.281) is equivalent to 

mB- 

As a first illustration of (I3.28L let p^’^^\p) = 2S{p), Then (I3.28P simplifies to a 
quadratic for G^, and from this use of the Sokhotski-Plemelj formula (11.81) implies 
the density of singular values is equal to 


1 / \ 1/2 

— ( — + (2a + 4)a)^ — ) , (3.29) 

where we have chosen f = 1/2 as in deriving (13.181) . supported on the region of the 
positive real axis such that the argument of the square root is positive. This result 
is well known; see e.g. [m Prop. 3.4.1]. 

As a second illustration, suppose p'^’(°)(^) = S{p — b) + S{p + h) so that the 
singular values of are all located at b. We then find that (|3.28l) gives a cubic 
equation for G^, 

g(z^(l - gf - a(l - g) - 9) = 1- g, (3.30) 

where we have set f = 1/2, and g := G^t jz. For a > 0, we can see from this that 
the density is supported away from the origin. The reasoning is that otherwise, for 
small z, g must behave like with 0 < a < 1, as would follow from its relation 

to G'^ and (I3.10|) . But this is incompatible with (13.301) unless a = 0. 


Equivalent viewpoints. The Green’s function (13.101) is the Stieltjes transform of 
the density of the singular values of the matrix W in (EB). The singular values of 
W also appear as the eigevanlues of the block matrix Y specified in (13.101) . Thus 
Y has n — m zero eigenvalues, m eigenvalues equal to the singular values of W, and 
m eigenvalues equal to minus the singular values of W; see e.g. (TT] Prop. 3.1.1]. 
Denoting the corresponding density, normalised to integrate to unitary, by p‘^^(x; r) 
we see that 

+ :^S(x). (3.31) 

2 + a 2 + a 

Hence, the corresponding Green’s function 


^ch/ 




r) := 


I 


p^'^ix] 


■ dx 


i+u-U 


z — X 
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is related to the Green’s function (jd.lOl) by 

(3.32) 

2 + a 2 + a z 

Substituting in p.l2l) gives the inhomogeneous Burger’s equation 

f ) + (2 + a)g^\z-, r)^g^\z-, r) + = 0 - (3-33) 

OT oz (2 + a) z'^ 

which is the form considered in [5] for the present setting. 

Another variant is to consider the eigenvalues of W^W, with W as in (13. ip . As 
already remarked below (IX^ . the eigenvalues are W^W are the squared singular 
values of W. Thus, with the corresponding density denoted by p^{X-,t), we have 

p^{X; t) = ^p-{Vx-, r). (3.34) 


The Green’s function for the eigenvalues of W is given in terms of p^ {X ; r) by 

P^{x;t) 


g^{z-T) = 


■ dx. 


1+ 


Recalling (I3.10p we thus have the relation 

5W(z^;f) = ^^. (3.35) 

Substituting in p.l2l) gives that g^ satisfies the Burger’s like equation 

T) + (a + 2zg'^{z- r)) ^^g^{z-, r) + g^{z; rf = 0. (3.36) 

This variant, modulo some rescaling, is the one given in US [7] in the present 
context. 

For completeness we present the solution to the initial value problem of both 
(I3.33P and (13.361) . which follow from (I3.28P by substitution. In the case of g^^, 
substituting (13.311) and (13.321) in (13.281) we find 


g^‘'(z;f) = 2Fch 


pch.(0)(^)^^ 


4“' - 


Fruaf 


1 


2z 


2 + 0, V z zFch — dr 


(3.37) 


where Fch = z — f(2 + d)g'^^. In the case of (13.361) . substituting p.35l) in p.28l) 
gives (I1.12p . 

With (11.121) of interest in mathematical statistics and thus as a stand alone result, 
let us show the method of characteristics discussed in Section [5] can be used to solve 
the partial differential equation (I3.36P directly. We first note (13.361) is an instance 
of (12.171) with replacements G = g^ and A = i, R = a + 2zg'^, C = —{g^Y. The 
DDEs (I2.18P and (12.191) read 

^f(a,/3) = 4 (3-38) 

^z(a,^) = a +2z(a,/3)g'^(a,/3), (3.39) 

■^g^{a,l3) = -g^(a,/3)^ 


( 3 . 40 ) 
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where (a,/3) are the transformed variables {z,f) and the abbreviated notation 
g^{a,j3) = {z{a, j3)-,T[a, j3)). The solution to these ODE’s with initial con¬ 

ditions z{a, 0) = a, f (a, 0) = 0 and {a, 0) = g^ (a) is 

= (3-41) 


z{a, /3) = 0^(1 -f g^ (a)/3) -k a(l -k g^ (a)/3)^, 

(3.42) 

1 + 5o («)/3 

(3.43) 

Since to proceed requires purely algebraic operations only, we suppress all of the 
arguments in what follows. We find a formula for gw,o from the last equation. 

„w „ 9 

9o i_gW^’ 

(3.44) 

and substitute it into (13.421) to obtain the formula for a. 


a = z{l-g^Pf-aPil-g^P). 

(3.45) 

JNext we turn to p.44p and hnd 


g^ = {l-g'^l3)g^ 

(3.46) 


which, including the formula (13.451) for a and (13.411) . gives an implicit solution of 

( 1 ^ 

g'^iz] f) = (1 - 2Tg'^{z; f)) g^ (^z (l - 2rg'^(z; r))^ - 2af (l - 2Tg'^{z; t))^ . 

(3.47) 


Since 


gW(z) = PV 



z — 


(3.48) 


with C i?'*' and is normalized to unity when integrating over g, > 0 

we see that (jl.l2p follows. 


4. Circular and Jacobi ensembles 

4.1. Circular ensembles. We now turn our attention to Smoluchowski-Fokker- 
Planck type dynamics of circular ensembles. These models arise when considering 
random (symmetric if /3 = 1, unrestricted if /3 = 2 or selfdual if /3 = 4) unitary 
matrices U of size N x N distributed according to the Haar measure; see e.g. m 
Ch. 2]. The Jacobian in this case reads 

J= (4.1) 

l<i<j<N 

so that eigenvalues lie on a unit circle. The diffusion is introduced based on 
the parametrization of U in terms of exponent of an Hermitan matrix. The joint 
eigenvalue PDF satisfies (12.51) with a drift term 

W = - log|e*'^^ 

and with the replacements {Ai} i-A- see e.g. [171 §11-2.1]. 


(4.2) 
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The log-gas system occupies the domain (j) € (—tt, tt] and is described by one and 
two body potentials 

V2[(j>,e)=-\ogW*-e^0\. (4.3) 

For the spectral density p° 

= (4.4) 

\i=i / 


in the global regime where r = Nt^ the hydrodynamical equation (HH) is equal to 


dp°{(t>;f) 

dr 


d(f> 




d(j) 


d(j)' log|e*-^-e*-^'|p°(0';T 


(4.5) 


We define a circular Hilbert transform as 

^o[/]((^):=ipV^#'cot(^^^)/(</>'), </>€/, (4.6) 

where / C (—7r,7r] and, since 9a;log|e“ — e®*'| = ^ cot (^^), the equation (14.51) is 
expressed as 


dp°{(l);f) 

dr 


A 

d(j) 


{p°{^-,T)no[p°{-,rm)- 


(4.7) 


To arrive at the final evolution equation, we introduce a circular Green’s function 


^“(^(r) = i^cot p°(y;r)dy, (4.8) 

which also satisfies (j2.6p and p.l5p with replacements G G° and p^ p°. Based 
on these properties, we again find the complex Burger’s equation 

^G°{z;t) + G°(z;T)^G°(z;f) = 0, (4.9) 

which is formally in the same form as the Gaussian case (12.161) . We can thus apply 
the same techniques to conclude that the solution of the initial value problem for 
this equation — initial spectral density — reads 

,,, ( 4 . 10 ) 

where is the initial support of The hydrodynamical equation (14.51) was 

first derived by Pandey and Shukla [301 Eq. (59)], using the hierarchy of equations 
satisfied by the dynamical correlation functions. The general solution (14.101) is given 
in [301 Eq. (63)]. Our main point here is therefore not a new result, but rather a 
common theme, namely the macroscopic hydrodynamical equation dHl). 

The particular case p°’^^'>{p) = 6{p) was studied in the context of QCD by [HIS]. 
Even though there is no closed form solution of (I4.10L several analytic features 
can be exhibited, including an effect analogous to that of the Gaussian ensemble 
evolution with initial condition (12.25^ : at a critical value to f two spectrum edges 
collide here being the left and right edges of the single interval of support. 
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4.2. Jacobi ensembles in trigonometric variables. We now move to the exam¬ 
ple of Jacobi ensembles. Consider a unitary (symmetric for (3 = 1, unconstrained 
by /3 = 2 or self dual for ^ = 4) matrix S of size (n + m) x (n + m) with n > m, 
divide it into 4 blocks 


5 = 


f '^nxn 

Xj'mxn 


nxm 


(4.11) 


and investigate singular values of sub-block t'. The corresponding Jacobian of this 
ensemble (see [HI §11.2.3]) reads 


n (4-12) 

j—1 


where a'= n — m+1 — ^ and \i G (0,1) denote non-zero singular values of t'. 

To obtain a Smoluchowski-Fokker-Planck equation for the joint PDF, we intro¬ 
duce new variables Xi = sin^, with (j)i G (0,7r). As was demonstrated in [T71 
§ 11 . 2 . 3 ], the new variables {4>i} permit an evolution ( 12 . 51 ) with a drift term 


/ m A y ^ A 

»' = -|y:iog»m“|-|y:iogcos>|- y: 

i—l i—1 l<j<k<m 


log 


Sin 


4'j ■ 2 

^ — sm — 
2 2 


(4.13) 


with b' = ^ and replacements {Ai} —>■ ^ m. When compared to the Jaco¬ 

bian (I4.12p . an extra b' term arise by transforming the measure dXi = ^ cos 
We read off one- and two body interactions from (I4.13L 


ViW 


a' . 2 b' 2 

logsm - --logcos -, 


and rewrite the latter 


= - In 


sin 



^2(0,0') 


log 




(4.14) 


so that the (j) —> —((> symmetry is evident. Accordingly, we form a spectral density 
of the form 

= + , (4.15) 

^ \i=l / 


normalized to unity when integrated over (j) G (0, tt) and even in </). Both the 
interaction term and spectral density has features present in the chiral spectral 
density (13.6p and the two body potential term ()3.5I1 . The current (lESD driving 
the time evolution of in the large m limit reads 




d f d 


d(t) \2 


■ log sin^ 


d(j)' log 


d 4> ^ 


= /(</'; o cot - -f - / dcj)’ cot 


sin — — sm — 
2 2 


r) . 




(4.16) 


where the d = ^ — 1, b' term has dropped out as subleading in the large m 
limit and the rescaled time parameter reads f = mr. The hydrodynamic equation 
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reads 

-^p\4’ir) = -^(p'^((?i;f)'Ho[a5(-) +p''(-; t)] (</>)), (4.17) 

where the Hilbert transform T-Lo was already defined in (14.61) . By using the proper¬ 
ties of Green’s function 

= (cot -foot )p°(y;f)d 2 /, (4.18) 


with h C (0,7r], we repeat the derivation of the complex Burgers equation (14.91) 
and obtain 


d 


dr 


(z; t) + hr cot - -f G (z; f) —G (z; f) 


d 


dz 


a G^{z; t) 
4 sin^ zl2 


(4.19) 


The equation has the same structure as the chiral Gaussian equation (1,4.1211 . and in 
fact reduces to that equation for small z. The underlying log-gas setup has therefore 
the same features — it consists of a fixed particle at <() = 0 of charge a and two 
mirror-like clouds for (j) G (—7r,0) and (j) G (0,7r) respectively. In the special case of 
vanishing charge a = 0, the resulting equation (14.191) coincides exactly with (14.91) 
obtained for the circular ensembles. 


5. Collective variables 

Collective variables is another approach to obtain the hydrodynamic equations 
(12.161) and (j4.12L The idea of collective variables was first introduced in plasma 
physics [9] and extensively applied to gauge theories [24]. Besides rederiving the 
aforementioned hydrodynamical equations, this method is suitable for obtaining 
asymptotic formulas for group integrals of Harish-Chandra/ltzykson-Zuber and 
Berezin-Karpelevich type. The former relate to the Gaussian ensembles whereas 
the latter appear in the chiral Gaussian ensembles. This relationship is the reason 
why we focus only on these two cases in this section, and don’t consider the circular 
or Jacobi spectral evolutions. 

5.1. Collective variables method. In the present context, one proceeds by trans¬ 
forming the Smoluchowski-Fokker-Planck equation (12.5p to new ’’collective” type 
variables A, 

Ai —>■ Aj({A}), {A} = (Ai,..., Aat) (5.1) 

where i = = l...iV'. These new degrees of freedom should a) use the 

symmetries of the system and b) have a well defined large JV limit. Typically 
iV' —>• oo from the beginning, and thus the particle system is treated as a fluid, so 
that the change is not bijective at least before taking the large JV limit. For the 
special case N' = N and N finite, the method corresponds to a bona fide variable 
change and was recently studied in the present context in [25] . The non-uniqueness 
of the collective variables means the aim is not an exact description in all regimes. 
However, since the new degrees of freedom conserve the symmetries, one expects 
to correctly reproduce certain macroscopic properties. 

Gonsider a general transformation X{q; {A}) with i index promoted to a variable q 
(i.e. N' —> oo) in a fluid approximation. This continuous case introduces functional 
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analysis by which the transformed Smoluchowski-Fokker-Planck equation 


5r7r^({A}) = L({A})7r^({A}), 


HW) 


h dX^ 


N 


E 


dW d 
dXi dXi ’ 


(5.2) 


obtained by writing Pt = exp>{—j5W)TTr in (12.51) . is transformed to a functional 
differential equation. 

According to eu, the function of {A} becomes a functional in the q variables 
7rT-({A}) = 7rr[X{q\ {A})]. Moreover, the Laplace-Beltrami operator L is re-expressed 
by the appropriate continuous chain rule 


a E 9 a.({a}) 9 d _ f ^A(g;{A}) 5 

9^3 ^ dXj d\ dXj J ® dXj SXiq)' 


where the X{q) are the new variables just as Xi in the discrete case. The transformed 
operator L = K + V reads 


_ 1 f a^A(g) 5 ^ f p p ^^( 9 ) i9A(p) S'^ 

dXf 5X{q)^^J ^ SA, dX, SX{p)6Xiq)' 

= _ fp _ 


(5.4) 

(5.5) 


where we supressed the {A} dependence in the coefficients. The transformed 
Smoluchowski-Fokker-Planck equation (15.2p is then 

5r7rT-[A] = (^.^[A] -I- F[A]^ 7fT-[A]. (5.6) 


5.2. Gaussian ensembles. In the case of Gaussian ensembles, the drift term W 
is given by (12:^ and the collective variable 


N 

A(g;{A}) = J^S(q-Xi) 
2=1 


(5.7) 


is the (non-averaged) one-point correlation function (see (12.71) 1. This choice is 
consistent with condition a) mentioned in the introduction to this section — it 
conserves the eigenvalue exchange symmetry. We calculate the kinetic part K with 
the help of the formula d\^X = —dqS{q — Xi), and the potential part V using 


1 



= PV f dp -^—6(X, 
Ji X,-p 


m), 


^ S{p - X,)S{q - Xj) = X{p)X{q) - S{p - q)X{p). 


We set an ansatz for the leading large N form of the joint PDF 


TTr = exp 




(5.8) 
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we find that the new functional St satisfies the evolution equation 


dr Sr = / dpX{p) 


- / dpX{p) 


1 52 . 

( dSr 1 

(3 dp"^ 1 

yJA(p) 

(n^ 

/ d 5 St 


\dpSX{p) 


5‘^St 


-nm) 


d 5 St \ 

dp 6X{p) I 


[dpSX(p)) ^^^^^^^dpSX(p)^ ’ 


(5.9) 


where TLlf] denotes the Hilbert transform (12.111) with supressed argument. 

We perform the large N limit of (15.91) by rescaling both the time Nt = t and the 
collective variable A = Np. In this limit, the first term on the RHS is subleading 
in N in comparison to the second and the time derivative on LHS. Ignoring this 
term, we obtain an equation for Sf in the Hamilton-Jacobi form. 


dfSf + 


dpp{p) 


1 / d SSf 

2 \dp6p{p) 


2 

-HpKp) 


d SSf 

dp 6p{p) 


= 0 , 


(5.10) 


where the position variable is p{p) and the conjugate momentum reads A(p) = 
This allows Sf to be interpreted as an action evaluated on a physical 
trajectory between p{p',T = 0) and p{p;t). The resulting Hamiltonian H = 
J dpp (i(cIpA)^ — 'H[p]dpA'^ needs a minor reformulation since it contains a prob¬ 
lematic Hilbert transform term. To this end, we invoke a canonical change of 
variables (p, A) —>• (p', A' = A -I- C) with C dependent only on p. This change 
leaves the Hamiltonian unaltered i.e. H'[p', A'] = H[p', A' — C\p']] and the action 
picks up a boundary term 


S'f = Sf + T\f - T| 


0 , 


(5.11) 


where the subscripts denote boundary terms evaluated at initial f = 0 and final 
time f. The generating function T is found to be 


T = -- 


1 


dpdqp{p)p{q) In \p - q\, 


6T 

5p{p) 


= C. 


The transformed Hamiltonian H' is 

1 


i/'[p',A'] = 


dpp'{p) {dpA\p)f - {n[p']{p)f 


(5.12) 


(5.13) 


for which the second term is reexpressed in term of p as 

J dpp'{p){'H[p']{p))^ = Y J dpp'{pf. (5.14) 

This identity is proved using the properties of the Hilbert transform 

/ mg] = - / gnn, rnmifw = inn? - 

valid for sufficiently well-behaved functions /, g [26] . From now on we drop the 
prime indices and the Hamiltonian (15.131) is finally 


H[P,^] = \J^PP - yP" 


( 5 . 15 ) 
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with corresponding action 

-S'f = ^ y di J dpp (^{dpAf + , (5.16) 

chosen so that S'f |f=o = 0 (otherwise Sf is unique only up to an additive constant). 
By the Hamilton equations dfP = ^,dfA = —^, the equations of motion read 

dfA + ^{dpAf = 

dfP + dp{pdpA) = 0. (5.17) 

Upon defining G± = ^iirp+dpA (cf. (12.141) 1. these formulas are exactly the complex 
Burger’s equation (12.161) . The construction goes similarly as before — equations 
(12.151) for G± are defined on the real line and induce a complex structure due to 
analytic properties of G. 

This is a well-known result of Matytsin [28], reproduced also by other authors 
[mdo]. Here we show how additionally the joint PDF function tt is asymptotically 
expressed in terms of an action related to the hydrodynamical system. 


5.3. Chiral Gaussian ensembles. For the chiral case we make suitable replace¬ 
ments (13.31) and W is defined in (13. 4p . The collective variable in this case is 

m m 

x{q, W) = ^ (5(q - Xi) + S{q -p Xi) = ^ 2|(7|d(g2 - xf), (5.18) 


a (non-averaged) one-point correlation function (see (13.61) 1. Our task is therefore 
to transform (lOD . Because the derivation is parallel to the Gaussian case, we give 
only some partial results. In calculating the transformed Laplace operator L, we 
use formulas 

|g| d 


^ 2|,|«(,> - X?)) = (2|,|i(,» - , 

S(q‘‘ - = PV f 


1 

Xq 


q dq 


-oo 111 - P 


2xi 


rr^ _ 


S{q^ - x]) = PV / 

J —C 


3 > 

dp,2\p\ 5(p2 - xl) 2 


kl - P 


5{q^-x^). 


We make a large m joint PDF ansatz 


= exp ( —’Arri'^S'^ 


a 


(5.19) 


which captures the rough degrees of freedom ~ rn? and trivial /3 dependence. The 
equation satisfied by S'^ reads 

- J ivm (jjp + 2jy) +«[.(,)](0)-y) + 


- / dpx{p) ( — 


f d \ 


d 5S% 
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Now we perform a m, n —>■ oo limit with njm, fixed. We set x = f = mr and 
find the first term on RHS subleading wrt. the second and the time derivative. The 
equation for S'? is again in the Hamilton-Jacobi form 


dfS%^ 


dpp%p) 



(d 6 Sl \ 
\dp 5 p’={p)) 


a' d 5 S% 
p dp Sp^{p) 


nip’^Kp) 


d SS^ \ 
dp 6 p^{p) j 


= 0 . 

(5.21) 


With the conjugate momentum A'^ = we transform {p^,A^) —>■ ((p'^)' = 

p^, = A‘^ -I- and the action is transformed accordingly 

(5|)' = 55+T|^.-rfo (5.22) 

with a generating function 

T‘= = -iy dpdq log \p - q\p^{p)p%q), ^ (5.23) 

Under this transformation, the new Hamiltonian reads 


i^[p^A=] 


1 

2 


dpp^ip) 


((d SS<f 
\ V(9^ Sp^{p) 


2 


- (hH(p))' 


2 d/ d SS^f 
p \ dp 5 p'^{p) 



(5.24) 


where a = ^ — 1 and we dropped the primed indices. As before, the term quadratic 
in Hilbert transforms is given in (15.141) whereas the linear term vanishes by two 
additional properties of % [55], 


n 


/(■) 


(a;) 


(•) 

^[/](0) = 0 


nn mm 

X X 

if f{x) = f[-x). 


(5.25) 

(5.26) 


Thus the final form of the Hamiltonian reads 

i^[p^ A^] = ly dpp^ [{d.A^^f - y(p=)^ + ySpA^ , (5.27) 

and the action is therefore 

= ll^diI dpp^ (^{dpA^? + y(p=)^ + jdpA^) ■ (5-28) 

The Hamilton equations of motion read 

5,A^ + i(apA^)2 + ^5pA^ = y (p^)^ 

dfP^ + dpip^dpA-^) + dp (^^p^ = 0. (5.29) 

We observe again how the chiral case reduces to Gaussian (12.161) when a = 0, and 
by defining G% = ^iirp^ + dpA^ we reclaim (13.121) by the arguments elucidated 
previously. 
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6. Asymptotic expansion of Harish-Chandra/Itzykson-Zuber and 
Berezin-Karpelevich integrals 

Collective variables were used by Matytsin to obtain large N expansion 
of the celebrated Harish-Chandra/ Itzykson-Zuber integral formula. Later works 
looked at the same task from both mathematical [2T] and physical point of view [10] . 
In this section we comment on this standard result and afterwards use analogous 
working to compute an expansion for the Berezin-Karpelevich type integrals [HlEo] 
arising in the chiral Gaussian ensembles. 


6.1. HCIZ-type integrals. We consider an integral 

= J{uUU)exp (^^Tr{UAU^B)^ , 


( 6 . 1 ) 


where matrices A, B are diagonal of size N x N and U are real orthogonal (/3 = 1) 
or complex unitary (/3 = 2). In the RMT context these integrals arise in con¬ 
nection with the Gaussian ensembles. For /3 = 2 an exact formula exists, found 
independently by Gharish-Handra [52] and Itzykson-Zuber [55] . 

To obtain large N asymptotic behaviour of (IQ) . we recall the definition dUD 
of joint PDF 


TTf 


j {UUU)Pf{X^°'>]ULU'^) 





^Tr(X(«))=j^ 


L A(0) 



The traces in this expression can be given in term of the collective variable (15.71) 
as TrL^ = ^ J d,pp^p{p;T) and Tr(A(°))^ = ^ J dpp^p{p-,T = 0) in accordance 
to their role as an initial and final densities respectively. On the other hand, the 
asymptotic form of tt,- was found in Section 15.21 as 

TTf ^exp(^-^A2(5f-T|^-rio)) (6.3) 

using (15.81) and (15.111) , and where we also added an arbitrary constant to the action 
5'f —)> S'f — 2T|o (recall the comment below (15.161) 1. The form of this constant is 
chosen so that Trf\f=o exp (/JA^Tjo). 

To arrive at an asymptotic expression for dSU), we fix the time r = 1 and rename 
the final p{p,T = 1) = Pfip) and initial p{p,T = 0) = Pi{p) densities 


1/3 (cr, a) 


^ CN,f=i exp 





i J dpdq(^pi{p)p,{q) + 


dp p'^{pi{p) + Pf{p)) + 

) ln|g-p| y 


(6.4) 


Now the main difficulty lies in finding a physical path joining initial Pi(p) and 
final pf{p) spectral densities and calculating the corresponding action Sf, which 
is specified by (I5.16|) . The former problem has been solved from our workings in 
Sections 1 and 2 for the initial condiitions (dU) and (|2.25|) . while the evaluation of 
Sf in the first of these is given in [10] . For a discussion of analyticity properties of 
(16.41) see [19]. 
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6.2. Berezin-Karpelevich type integrals. We now turn to the asymptotic for¬ 
mula for an integral of Berezin-Karpelevich type defined as 

J 0 {A, B)= J{UUU){VUV) exp (^^Tr {VA'^U'^B + B^UAV"^)^ , (6.5) 

where A,B are nx m diagonal matrices and U,V are real orthogonal (/3 = 1) or 
complex unitary (/3 = 2) matrices of sizes nxn and mxm respectively with n > m. 
These integrals arise in studying chiral/Wishart/Laguerre type ensembles. 

In the (3 = 2 case, an exact formula was rediscovered in [20] and originally 
calculated by Berezin and Karpelevich [3]. To obtain an asymptotic expression for 
(3 = 1, 2 we recall the chiral joint PDF (13.2p . 


: J {UUU){VUV)Pf{Z^°'>;ULV^) 


C 




Jp 


N- 


L \ 


( 6 . 6 ) 


We introduce the normalized collective variables to the Gaussian terms TrL^L = 
^Jdq q^p^{q]T) and Tr(Z(°) )tz(0) = ^ J dq q^p%q;T = 0). An asymptotic form 
of LHS was found in Section 15.31 as 



along with adding a constant — 2T^^. By comparing (16.61) and (16.71) . for 

fixed time r = 1 we have an asymptotic formula 


Jp (p/, P/) C'A,f=i 6xp +2 J (Pi(p) + P/(p)) + 

+ l f dpdq {p’({p)p'(iq) + p^p)p}iq))ln\q - p\ 


( 6 . 8 ) 


where initial and final densities are denoted as p‘^(p;T = 0) = Pi{p) and p‘^(p-,T = 
1) = p^(p) respectively. As in the case of (16.4p . we comment that to obtain the 
asymptotic formula for prescribed initial p\ and final p’j densities, it is necessary 
to evaluate the action S’; on a physical trajectory connecting these two spectral 
densities. We note that such trajectories are given for particular initial conditions 
below p.28p . 
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